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Alpha decays in the EXO-200 detector are used to measure the fraction of charged ®®®Po and 
®®®Bi daughters created from alpha and beta decays, respectively. ®®®Rn alpha decays in liquid 
xenon (LXe) are fonnd to produce ®®®Po"'" ions 50.3 ± 3.0% of the time, while the remainder of 
the ®®®Po atoms are neutral. The fraction of ®®®Bi+ from ®®®Pb beta decays in LXe is fonnd to be 
76.4±5.7%, inferred from the relative rates of ®®®Po and ®®®Po alpha decays in the LXe. The average 
velocity of ®®®Po ions is observed to decrease for longer drift times. Initially the ions have a mobility 
of 0.390 ±0.006 cm®/(kV s), and at long drift times the mobility is 0.219 ±0.004 cm®/(kV s). Time 
constants associated with the change in mobility during drift of the ®®®Po^ ions are fonnd to be 
proportional to the electron lifetime in the LXe. 
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I. INTRODUCTION 

EXO-200 discovered the two-neutrino double beta de¬ 
cay in ®®®Xe [I] , and subsequently made a precision mea¬ 
surement of the half-life of the decay [2]. The experi¬ 
ment is currently searching for neutrinoless double beta 
decay {OiyPjL), and has produced limits on the half-life of 
this process in excess of 10®® years [ail] using 110 kg of 
liquid Xe (LXe) enriched to 80.6% in the isotope ®®®Xe 
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Figure 1. (Color online) Partial diagram of the uranium se¬ 
ries. The alpha particle energies are shown along with their 
corresponding relative intensities. The beta decays are shown 
with their branching ratios. 


[3. Analyses of alpha decays were used to establish the 
contributions of internal backgrounds (e.g. from the ura¬ 
nium series) to the QvfijS measurement [6j. In this pa¬ 
per, these alpha decays are used to make novel mea¬ 
surements of ion fractions resulting from alpha and beta 
decay and to determine the ^^®Po ion mobility in LXe. 
These measurements may also be of interest to other ex¬ 
periments that utilize a noble liquid TPC, such as dark 
matter searches and long-baseline neutrino oscillation ex¬ 
periments. Radon and its by-products are also found in 
these detectors. 

The relevant portion of the uranium series is shown 
in Fig. [1] This part of the series contains three alpha 
and two beta decays that occur within a relatively short 
time period. Using the EXO-200 detector it is possi¬ 
ble to match many ^^^Rn decays to their daughter ^^®Po 
decays. These matched pairs are then used to measure 
the ion fraction and ion mobility of ^^®Po resulting from 
the alpha decay. Finally the activities of the three alpha 
decays,^^^Rn, ^^®Po, and ^^"^Po, are used to determine 
the ion fraction of ^^^Bi produced from beta decay. 

The fraction of ions created from beta decays and the 
behavior of ions in the detector is of particular interest 
for nEXO, the next-generation EXO experiment, which 
is expected to contain 5 tons of LXe enriched in ^^®Xe. A 
possible second phase of nEXO may include a system to 
“tag” the ^^®Ba daughter of the Ov/SP decay in order to 
eliminate all background except a negligible contribution 
from the two-neutrino mode 13 [5]. There are several 
important unanswered questions relevant to the imple¬ 
mentation of Ba tagging: 1) whether the ^^®Ba daughter 
is mainly an ion or a neutral atom, 2) what the charge 
of the ion is in LXe, and 3) how the Ba will move in the 
detector over time scales required to capture it. One way 
to probe the behavior of ions and atoms in the detector is 
to track events as they decay through the uranium series. 

To our knowledge there are no previous determina¬ 
tions of the daughter ion fractions of alpha or beta decay 
within liquid noble gases. In related work with super- 



Figure 2. Cutaway view of the EXO-200 detector highlighting 
several important components of the TPC. 


fluid helium, the ion fraction of ^^®Rn recoils into the 
liquid from ^^^Ra alpha decay occurring on a surface has 
been measured to be 5.36 ± 0.13% and 1.04 ± 0.06% for 
two different surface geometries [S] , and the surviving ion 
fractions for high energy external beams of ^^B, and 
®Li have been reported as 30%, 10%, and 35%, respec¬ 
tively m- The ^^®Po ion fraction from ^^^Rn decays in 
helium and argon gas has been measured to be 59 ± 4% 
and 48 ± 5% respectively at atmospheric pressure m- 
The ion fraction of ^^®Po has also been measured in air 
at atmospheric pressure to be 87.3 ± 1.6% m, and 88% 

m- 

Mobility measurements have been reported for four al¬ 
kaline earth ions in LXe including Ba’*' as well as 
Tl'*’ [m [15] and Th^ [16] . These mobilities range from 
0.133 cm^/(kV s) to 0.280 cm^/(kV s). 


II. EXO-200 DETECTOR 

EXO-200 is described in detail elsewhere |S]. To 
summarize, the detector (shown in Fig. consists of 
two symmetric back-to-back cylindrical time projection 
chambers (TPC) with a common mesh cathode. The 
TPCs are labeled TPCl and TPC2 for the +z and —z 
sides of the cathode respectively, (where 2: is along the 
axis of symmetry). The detector is filled with LXe that 
is isotopically enriched to 80.6% in ^^®Xe, with the re¬ 
mainder being primarily ^^^Xe. The density of this LXe 
is 3.03 g/cm^. LXe serves as both the double beta de¬ 
cay source and the detection medium, converting energy 
deposits into ionization and scintillation. The LXe tern- 
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perature in the detector is 166.6 ± 0.2 K, and nominal 
pressure is 146 ± 6 kPa. 

The central cathode is set to a potential of —8.0 kV. 
At each anode two crossed wire planes are used to mea¬ 
sure the drifted ionization. The potential of each plane 
is set such that the charge drifts past the first wire plane 
(V-wires = —780 V), inducing a signal, and is collected 
on the second wire plane (U-wires = 0 V). These ion¬ 
ization signals are used to determine the event’s x and 
y position(s). There are 38 U- and 38 V-wire channels 
in each TPC, each channel consisting of three adjacent 
wires connected together. The wire-to-wire spacing in 
both wire planes is 0.3 cm. The total drift distance from 
the cathode to the V-wire plane is 19.2 cm. There is a 
0.6 cm gap between the V- and U-wire planes. Behind 
the U-wires, by 0.6 cm, is a plane of large area avalanche 
photo diodes (APDs) which detect the scintillation light 
|17j . These APD planes are biased with —1400 V. A 3D 
simulation of the TPC has determined the drift held in 
the analysis region to be 380 V/cm, with some position 
dependence (±5 V/cm). 

The data acquisition system (DAQ) continually sam¬ 
ples signals from each U, V and APD channel at 1 MHz. 
When a U-wire, APD, or APD-sum goes above appro¬ 
priate thresholds, the DAQ records a 2048 ys long wave¬ 
form, with 1024 pre-trigger samples. Later, during data 
processing, all waveforms are scanned for signals. APD 
and wire signals are then grouped together using spatial 
and temporal information. The time difference between 
scintillation and ionization signals is converted to the 2 
position of the event, using the electron drift speed of 
1.71 mm//is [2]. If two events occur within the maxi¬ 
mum electron drift time (116 /is), a clustering algorithm 
associates the charge signals to the scintillation signals 
that occurred closest in time. Further information about 
data handling and reconstruction of events in EXO-200 
can be found in [2]. 

High purity conditions are maintained in the detector 
to maximize collection of charge and light. This purity 
is achieved by continually recirculating xenon gas in a 
purification system consisting of a heater, a pump m, 
getter purifiers, and a condenser. The xenon flow is in 
the -2 direction with respect to the detector coordinates, 
with the input near the TPCl anode, and output near 
the TPC2 anode. The purity of the LXe in each TPC is 
monitored using a calibration source deployed every few 
days. Purity is measured in terms of the electron lifetime, 
as determined from the attenuation of the charge signals 
produced by the calibration source as a function of drift 
distance. The electron lifetime is typically maintained 
between 3 and 5 ms. A xenon feed and bleed system is 
used to fill and empty the detector, as well as to make 
adjustments to the xenon system pressure. Occasional 
pauses in recirculation, for various operational reasons, 
cause temporary drops in the LXe purity. These recir¬ 
culation pauses often coincide with feed events, in which 
stored Xe is introduced into the system, which further re¬ 
duces the electron lifetime since this Xe has not had the 


benefit of constant purification. Feed events also intro¬ 
duce ^^^Rn into the detector from the gas storage system. 
The source of this ^^^Rn is unknown, but is likely due 
to emanating component (s) or from known small leaks 
to atmosphere (or both) in the gas storage system. The 
xenon gas purity in EXO-200 has been measured for some 
contaminants in the recirculation system (gas phase) us¬ 
ing cold trap mass spectrometry during operation m, 
however these measurements are not necessarily repre¬ 
sentative of the contaminants in the liquid phase, nor 
are they exhaustive. 

The experiment is installed underground at the Waste 
Isolation Pilot Plant near Carlsbad, NM. The site has an 
overburden of 1585 meters water equivalent which par¬ 
tially shields the detector from cosmic rays. An active 
muon veto system is installed around the detector to help 
reject cosmogenic backgrounds. 

III. ANALYSIS 

This analysis was performed using the data set from 
the latest OvPP result [4] with the addition of runs near 
xenon feed events. These additional runs were not used in 
d] due to the rapidly changing electron lifetime. Never¬ 
theless, they are useful here because they have a higher 
^^^Rn content, and this analysis, relying primarily on 
scintillation signals, does not require high accuracy for 
charge energy. Time periods following muons detected 
by the veto system or the TPC are not removed since 
cosmogenic backgrounds are unlikely to be confused as al¬ 
pha decays. However, events containing muon tracks and 
events identified as noise are removed from the data set. 
The total live-time used for this analysis is 572.8 days. 

Several corrections are applied to the data including 
wire-channel gain, electron lifetime, and V-wire shield¬ 
ing inefficiency (see [5] for more details). Eor events that 
are fully reconstructed (have a U-, V-wire, and a scintil¬ 
lation signal), the APD gain variations and the position 
dependence of the scintillation light collection are cor¬ 
rected using a light map derived from calibration source 
data. Events missing V-wire signals cannot be fully re¬ 
constructed in x-y position, and therefore cannot be light 
map corrected. All of the events used in this analysis have 
at least a scintillation and U-wire (charge collection) sig¬ 
nal, which is sufficient to determine the z position of the 
event. 


A. Identifying alpha decays 

Alphas interacting in the LXe are readily distinguished 
from beta and gamma interactions due to the large 
amount of scintillation relative to ionization. This differ¬ 
ence is caused by the higher ionization density for tracks 
created by alpha particles, resulting in more recombina¬ 
tion. This is demonstrated in Fig.j^where — )■ ^^^Po 

(Bi-Po) coincidence events are plotted as scintillation 
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Figure 3. (Color online) Raw scintillation versus electron 
equivalent charge energy (keVee) of Bi-Po events. Different 
colors distinguish between Bi and Po on surfaces and in the 
bulk LXe. Surface events are defined to be within 10 mm of 
the cathode or anodes. The inset shows a histogram of the 
time delay between Bi-Po pairs. 


counts versus charge energy. The decay (emitting 

a beta and several gammas) and (alpha decay) are 

easy to identify since the events occur in rapid succession 
(Ti /2 = 163.8 fis), usually in the same DAQ frame. Here 
“bulk” coincidences (blue and black) are identified by a 
^^^Po alpha decay with a z position between 10 mm and 
182 mm from the cathode, while “surface” coincidences 
(magenta and green) are within 10 mm of either the cath¬ 
ode or V-wire surfaces. Alpha decays on or near the 
cathode or anodes are observed to have a higher charge 
energy, and a correspondingly lower scintillation energy, 
than those which decay in the bulk xenon due to the 
higher local electric field. 

A scintillation versus charge histogram of the whole 
data set, highlighting the regions of interest with alpha 
events, is shown in Fig. for both raw and light map 
corrected scintillation counts. In the fully reconstructed 
and light map corrected data (upper right), the posi¬ 
tion of alpha events in the bulk (region B) and near 
surfaces (region C) is clearly separated, which is shown 
in the lower plots. The three alpha decays in the ura¬ 
nium series are identified with as three separate peaks in 
Fig.|4] region B, with mean scintillation of 38766, 42377, 
and 54254 counts for the ^^^Rn, ^^®Po, and ^^“^Po de¬ 
cays respectively. These means are found to be linear to 
within < 0.03% with respect to the alpha energies. There 
are significantly fewer alpha events in the light map cor¬ 
rected data, due to missing V-wire signals. To maintain 
full statistics, these events are analyzed using the un¬ 
corrected scintillation and U-wire signals, removing the 
requirement for a V-wire signal. 

The rate of alpha decays of ^^^Rn and ^^®Po from the 
region A highlighted in Fig. (upper left), as a func¬ 
tion of time spanning the data set, is shown in Fig. 
along with the average electron lifetime and the xenon 
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Figure 4. (Color online) Upper-Left: raw scintillation counts 
versus charge energy histogram. The highlighted region (a) 
is used to monitor activity of the ^^^Rn chain in the LXe. 
Upper-Right: fully reconstructed and light map corrected 
data. Lower-Left: positions of events in region (b). Lower- 
Right: positions of events in regions (c) (dark blue) and (d) 
(red). 

recirculation flow. The rate of ^^^Rn and ^^®Po alpha 
decays points to two sources: the constant emanation 
from some component(s) in the recirculation loop, and 
injection from the high pressure xenon gas supply which 
is introduced during xenon feed events. 


B. Counting alpha decays 

While fully reconstructed alpha events have good scin¬ 
tillation energy resolution, the V-wire detection efficiency 
is low (< 30%) and not well characterized at such small 
charge energies. Thus, for identifying alpha decays in 
the bulk LXe, only the scintillation and U-wire signals 
are used. The detection efficiency is > 99.95% for U-wire 
and APD channels for ^^^Rn decays (the lowest energy 
alpha decay considered). 

Alpha like events are selected if the event contains at 
least one scintillation signal between 25k and 66k counts. 
Then the event then must pass one of three criteria: (a) 
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Figure 5. (Color online) Top: rate of alpha decays in the re¬ 
gion defined by Fig. (A), primarily from ^^^Rn and ^^®Po 
decays. Middle: average electron lifetime measured in the 
TPCs. Bottom: xenon recirculation flow rate. Red lines indi¬ 
cate when xenon feed events occurred, wider lines correspond 
to larger feeds. 


the event has a single charge cluster and a single scintil¬ 
lation cluster (i.e. ^^^Rn or ^^®Po), (b) the event has two 
scintillation clusters (Bi-Po), or (c) the event has one or 
more charge cluster(s) collected from the Bi decay before 
the only scintillation cluster (Po) (i.e. Bi scintillation 
is below threshold). For Bi-Po events (b and c) a time 
window cut is applied which only selects events that are 
within 120 and 800 fis apart. This ensures that all charge 
clusters from the first decay (Bi) have drifted out of the 
TPC before the second decay (Po) occurs, and that the 
charge cluster of the second decay is collected prior to 
the end of the waveform. The efficiency of counting Bi- 
Po events in this time window is 55.7 ± 7.2%, which is 
the fraction of exponential decay in the time window with 
a small correction due to events of type (c) which have 
an uncertainty due to lacking the Bi scintillation signal. 
This factor was determined by applying the event selec¬ 
tion criteria to Monte Carlo events with realistic thresh¬ 
olds and noise. 

In Fig. (top), the raw scintillation counts of alpha 
events found in this manner are plotted versus the z po¬ 
sition. It is seen that the alpha decays of ^^^Rn, ^^®Po, 
and ^^"^Po are partially resolved. The resolution of the 
raw scintillation is degraded towards the anodes and the 
cathode, where x-y position dependent effects are more 
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Figure 6. (Color online) Top: raw scintillation signal z po¬ 
sition dependence for ^^^Rn, ^^®Po, and ^^‘*Po alpha decays. 
Vertical lines designate range of the data used (30 < |z| < 
150 mm). Bottom: histogram of the ^-corrected scintillation 
signals in TP Cl with fits to each decay. 


Table I. Number of decays determined by the fits of the scin¬ 
tillation signals, per TPC, in the range 30 < \z\ < 150 mm. 
Note the efficiency for counting ^^^Po is 55.7 ± 7.2%. 


Decay 

TPCl 

TPC2 


8036 ± 189 

8069 ± 188 

218 po 

5549 ± 106 

5286 ± no 

2i4po 

518 ±23 

521 ± 23 


significant. A linear z-dependent correction is applied 
to events in the range 30 < |z| < 150 mm in each TPC. 
The resulting histogram of these events is shown in Fig.[^ 
(bottom) for TPCl (with a similar histogram for TPC2). 
These peaks are then fit with three left-skewed Gaus¬ 
sian functions with a common skew factor and widths 
constrained by the ratio of their median energies. The 
resulting fits are integrated to determine the number of 
events associated with each decay in the analysis region 
(shown in Table |l]). The parameter errors reported by the 
fit are used to establish the uncertainty on the number 
of decays of each alpha. The left skewed Gaussian func¬ 
tion was chosen based on the best overall after also 
considering Gaussian and Grystal Ball functions. 

The number of bulk events associated with subsequent 
decays in the series is seen to decrease. This is because 
each decay in the chain has a probability of producing a 
positively ionized daughter that drifts towards the cath¬ 
ode and out of the analysis volume. The fraction of ions 
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Figure 7. (Color online) Raw scintillation counts versus a 
position of coincident pairs assigned as ^^^Rn (red) and ^^®Po 
(black) in the region 20 < \z\ < 172 mm. 


produced by each decay is discussed in Section |IV C| and 
Section IIVDI 


C. coincidences 

Using both the spatial and temporal information pro¬ 
vided by the TPC, a ^^®Po decay can often be associated 
with its parent ^^^Rn decay. Such delayed coincidences 
provide information about the ^^®Po charge state and the 
motion of individual daughter atoms or ions in LXe. 

Events used to search for ^^^Rn-^^®Po (Rn-Po) coinci¬ 
dences have a single charge cluster, a raw scintillation 
signal between 28k and 50k counts, and are between 
20 < 1^1 < 172 mm. The pairing algorithm requires 
that events occur within 3 minutes of each other and have 
charge collected on the same or on an adjacent U-channel, 
allowing for small transverse motion due to diffusion or 
xenon motion. In the case of two fully reconstructed 
events, they are also required to be detected by the same 
or an adjacent V-wire channel. Alpha candidate events 
that have more than one possible coincidence match are 
removed from the analysis (< 1% of the total). Coin¬ 
cidence events are assigned as ^^^Rn and ^^®Po decays 
based on their order in time. This technique is possible 
because of the very low rate of alpha decays in the detec¬ 
tor. Fig. shows the raw scintillation signal versus the 
z position of all paired Rn and Po events. The energy 
differences between Rn and Po populations corroborate 
the decay assignment based on time ordering. The total 
number of Rn-Po pairs found in the data set is 6507, with 
277 of these having both decays fully reconstructed and 
2083 with one decay fully reconstructed. 

The EXO-200 Geant4-based Monte Carlo [5] was used 
to determine the probability that two alpha decays ran¬ 
domly placed in the LXe would have their charge col¬ 
lected on the same or adjacent U-wire and therefore be 
found in coincidence using the pairing method described. 



Figure 8. (Color online) Scatter plot of ^^®Po drift distance 
versus time between the ^^^Rn and ^'^®Po decays. Displace¬ 
ment (Az) is defined as positive when movement is towards 
the cathode. 


This simulation has a detailed model of the detector, in¬ 
cluding the wire channel layout, and determined that the 
spatial requirement rejects 95.9% of random (false) coin¬ 
cidences. 

The rate of false coincidences can also be estimated 
from the data, using the fully reconstructed candidate 
events. For this data set, matching done with only the U- 
channel information finds an additional 16 coincidences. 
Further investigation of those 16 pairs identifies 3 with 
separations in the x-y plane of > 100 mm, while the 
remaining 13 events have light map corrected scintillation 
energies and separations in the x-y plane (< 30 mm) that 
are compatible with valid coincidences. The number of 
false coincidences is therefore estimated at < 2% for the 
Rn-Po pairs matched using U-wires only. 

The position difference (Az) versus the time difference 
(At) for these Rn-Po coincidences is shown in Fig. 
where Az is defined as positive for a displacement to¬ 
wards the cathode. Two primary populations are ob¬ 
served which are attributed to ^^^Rn that: (a) produced 
neutral ^^®Po atoms which only moved slightly, and (b) 
produced positive ^^®Po ions which drifted towards the 
cathode. A histogram of the mean drift velocity of these 
^^®Po ions and atoms is shown in Fig. The mean ve¬ 
locity of ^^®Po ions is significantly broader than for the 
^^®Po neutrals. The cause for this broadening of ions is 
discussed in section ITV Bl 

Matching ^^®Po and ^^"^Po decays has also been at¬ 
tempted with the method described here. Due to the 
longer decay time, accurate assignment is difficult with- 
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Figure 9. (Color online) Histogram of the mean velocity 
of ^^®Po ions and atoms extracted from ^^^Rn-^^®Po coin¬ 
cidences. 


out improvements to the V-wire detection threshold 
and/or to the matching algorithm. 


IV. RESULTS AND DISCUSSION 

A. Liquid xenon flow 


The flow of LXe in the detector can be measured us¬ 
ing the Rn-Po coincidence events that produce neutral 
^^®Po atoms. The detector is divided into 12 equal vol¬ 
umes, with 6 bins in 2 and 2 bins in the transverse di¬ 
rection using the U-wire position. Velocity histograms of 
the neutral pairs in each volume are fit with a Gaussian 
to determine the mean velocity and error. The results, 
shown in Fig. 10 indicate non-uniform flow within the 
detector. The observed mean velocity of ^^®Po atoms in 
the active volume of TPCl is 1.9 ± 2.6 x 10“^ mm/s, 
while the observed mean velocity in the active volume of 
TPC2 is —26.4 ± 2.3 x 10“^ mm/s (green bands). Note 
that, because of our convention for Az and v in Figs. 
and 1^ these numbers imply motion in the same direc¬ 
tion. For comparison, a velocity of 3.6 x 10“^ mm/s 
in the —z direction would be obtained for a recircula¬ 
tion rate of 14 SLPM (gas), assuming that the flow was 
uniform across the « 20 cm detector radius (red line). 
The average motion is aligned with the direction of LXe 
recirculation, but it is quite non-uniform, implying sub¬ 
stantial flow outside of the active regions of the TPCs 
(i.e. outside the teflon reflector or hexagonal wire plane 
regions). Note that the average flow direction is the same 
in the bottom of both TPCs, but it is in opposing direc¬ 
tions in the top regions of the two TPCs. In either case 
fluid flow is significantly less than the ion velocity (of 
order 1 mm/s). 



Figure 10. (Color online) Top sketch: the “upper” and 
“lower” volumes in each TPC are defined using the U-wires 
as shown in the hexagonal active volume. Black circle shows 
the internal diameter of the detector, and red circle shows 
the position of the teflon reflector. Green arrows show where 
the LXe flows into and out of the detector at the end caps, 
near the anodes. Bottom: the mean velocity of neutral ^^®Po, 
measured in the different volumes of the EXO-200 detector. 
The green band shows the average velocity of all events in 
each TPC, and the red line shows the expected velocity due 
to the LXe recirculation. Note that the vertical axis on TPCl 
is reversed to account for the convention on the sign of the 
velocity in which velocities are defined as positive for motion 
towards the cathode. 


B. 


^Po ion drift velocity and mobility 


The broadening of the velocity distribution of positive 
^^®Po ions relative to that of neutral ^^®Po atoms in Fig.[^ 
can be understood by plotting the drift time (At) versus 
velocity, as shown in Fig. El It is seen that the average 
ion velocity decreases with At. A possible explanation 
for this effect is that the ^^®Po ion initially moves with 
a higher drift velocity Ui, and while drifting a reaction 
or charge transfer occurs resulting in a larger molecular 
ion or a reduction of the charge of the ion. After this 
reaction the ion moves with a lower drift velocity, V 2 - 
Since the velocity is measured using the initial and final 
decay positions and times, the observed velocity is the 
time average spent at vi and V 2 - 

It is interesting to compare the velocity distribution for 
data sets with different ranges of electron lifetime (re). In 
Fig. 12 it is observed that the ion velocity distributions 
shift towards lower values as Te decreases, suggesting that 
the presence of impurities affects the ion transport. A 
Monte Carlo model with five input parameters (ui, ^ 2 , 
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Figure 11. (Color online) The velocity versus time between 
Rn-Po coincident events (positive is towards the cathode). 
The TPC in which the events occur is indicated by the color 
of the points. 



Velocity (mm/s) 


C, N, D) is used to describe the ion velocity data shown 
in Fig. 12 In addition to vi and V 2 (described previ¬ 
ously), the parameters C and N are the ratios of the ion 
reaction and neutralization time constants to and D is 
an effective 1-dimensional diffusion constant. To match 
the actual Te distribution in the data, the Te in each simu¬ 
lation is randomly selected from the data set of Te values 
of all Rn-Po coincidences. The initial charge state of the 
^^®Po is randomly chosen with probabilities that match 
the observed ratio of ions to atoms (see Section IV CI. 
The simulation then generates a random decay time (At), 
reaction time, and neutralization time from their respec¬ 
tive exponential distributions, and the initial 2 position 
from a uniform distribution. The Az of the ion drift 
is then calculated. Diffusion and an additional term, 
5z, representing the 1 mm uncertainty of z position (for 
details see m) is added to the final position as a nor¬ 
mally distributed offset with variance 2D ■ At + Sz^. The 
Monte Carlo generates 20-million simulations, which are 
histogrammed and scaled by the ratio of coincidences to 
simulations. The model parameters are found by mini¬ 
mizing the combined of all four histograms, and un¬ 
certainties are found by manually profiling each variable 
at this minimum. The minimum has a of 73.9 with 
63 degrees of freedom. The best fit parameters values 
with 1-sigma errors are vi = 1.48 ± 0.01 mm/s, V 2 = 
0.83 ± 0.01 mm/s, C = 12600 ± 660, N = 6.0lf;? x 10®, 
D = 0.61 ± 0.04 mmVs. 

The value of D is two orders of magnitude larger than 
the expected diffusion coefficient of an ion in LXe, so 


Figure 12. (Color online) Histograms of the average velocity 
of positive ions separated into four ranges of electron lifetimes. 
A model (solid lines) is used to fit the data (points with error 
bars). The dashed lines indicate the velocities ui and V 2 of 
the fit. 


this parameter presumably represents the effect of xenon 
motion rather than proper diffusion. The agreement be¬ 
tween this model and data supports the reaction model 
with a reaction time constant proportional to Tg. Another 
scenario to consider is that the initial species is 2 i 8 pQ-i-+ 
and the reaction is to ^^®Po^. This charge transfer could 
be aided by an impurity of lesser ionization potential than 
2i8po-i-^ However, this would require the survival of only 
charge 0 and 2-|- daughters after the decay, with no sig¬ 
nificant fraction of 1-|- daughters, an a-priori unlikely sce¬ 
nario. The effect of charge-induced fluid motion by the 
drifting xenon holes has also been considered, but esti¬ 
mates indicate that it is negligible. Using the drift elec¬ 
tric field in the detector, 380 ± 5 V/cm, the mobilities of 
the two drifting species are 0.390±0.006 cm^/(kV s) and 
0.219±0.004 cm^/(kV s), for vi and V 2 , respectively. The 
slower mobility {V 2 ) is similar to mobility measurements 
of other atomic ions, while the initial mobility (ui) is 
significantly higher than any ion mobility measurements 
found in the literature (with the exception of xenon holes, 
which are understood to be Xe^’*’ ions with charge motion 
predominantly by resonant charge transfer |20j). 
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Electron lifetime (ms) 


Figure 13. (Color online) Left: velocity histogram after apply¬ 
ing the At cut on all of the Rn-Po coincidences. The shaded 
region (velocity > 0.5 mm/s) is integrated to determine the 
number of ^^®Po ions. Right: the ion fraction found for dif¬ 
ferent electron lifetime ranges. 


C. Alpha decay: ion fraction 


component for the decays in bulk LXe originates from 
the ^^®Po decays that result in neutral atoms. This is 
expressed as 1 — In addition, the fraction of ionic 
that decay while drifting, epb, gives a contribution 
fa ■ epb to the measured number of ^^^Pb decays. The 
ion fraction, /^, from the alpha decay of ^^®Po, may be 
different from /„, the ion fraction from alpha decays of 
^^^Rn. However, it is reasonable to approximate with 
the value of fa measured in Section IV C| since the decay 
energies are similar and the ionization potentials (8.41 eV 
for Po and 7.41 eV for Pb) only differ by 1 eV. The factor 
epb = 0.088 ± 0.023 is calculated using the Monte Carlo 
model and parameters in section IV B[ with the half life 
of ^^®Po and an initial z distribution of ^^®Po calculated 
from a previous simulation of Rn-Po drift, which started 
with a uniform z distribution of ^^^Rn from the anode 
to the cathode. The number of ^^"‘Po decays occurring 
in the bulk LXe is similarly related to the number of 
^^"‘Pb events, with as the ion fraction of ^^^Bi and ion 
contribution factor epi = 0.108 ± 0.028 calculated with 
another iteration of simulations of ^^"‘Pb. The ratio of 
^^"^Po events to ^^®Po events in the bulk, 

is therefore equal to 


Using the Rn-Po coincidences it is possible to calculate 
the fraction of ^^^Rn alpha decays that produce ionized 
^^®Po (/a). To have an equal probability of selecting neu¬ 
tral and ion daughters, a cut on At is applied to ensure 
that no ions can drift outside the analysis volume. This 
cut is At < |zRn|/Pm, where zpn is the position of the 
Rn decay and Vm is a conservative maximum velocity of 
2.5 mm/s. Fig. 13 (left) shows the velocity distribution of 
all coincidences passing this unbiased cut. The ion frac¬ 
tion is the number of ions (entries with velocity between 
0.5 mm/s and 2.5 mm/s) divided by the total number of 
entries. As shown in Fig. [I^ (right), /„ is independent of 
the electron lifetime within statistical uncertainties. It is 
found that fa = 50.3 ± 3.0%, where the error is statisti¬ 
cal, and the uncertainty on the separation of the neutral 
and ion populations is negligible. To our knowledge this 
is the first measurement of the daughter ion fraction of 
an alpha decay in a noble liquid. 


D. Beta decay: ion fraction 

The ion fraction of ^^^Bi (fp), the daughter of ^^^Pb 
beta decay, can be estimated using the observed activity 
of the ^^®Po and ^^"^Po decays in the analysis volume (see 
Table 1^. As seen in Fig[^ the chain following ^^®Po is 

218po i!, 214p^^ 4 . 214 pq) 210pb. 

Here the ^^^Bi and ^^^Po decays are lumped together 
because negligible ion drift occurs during the short ^^^Po 
half life. 

Since the ^^^Pb ions are swept to the cathode in a 
time that is short compared with their half life, the main 


= (1 - /a + /QePb)(l - //3 + fpiBi)- 

From Table|^and the ^^^Po efficiency factor, = 

0.172 ± 0.023. In addition to the statistical uncertainties 
in Table the difference of the measured ratios in the 
two TPCs and the mean was included as a gauge of sys¬ 
tematic errors in the skewed Gaussian fits. Under the 
assumption that = fa, the equation above can be 
solved for fp, finding fp = 76.4 ± 5.7%. 

The same prescription can be followed to relate the ac¬ 
tivity of ^^^Po to that of ^^^Rn, with the addition of one 
more step (^^^Rn-^^®Po). The motivation for this is to 
check that the left-skewed Gaussian fit does not system¬ 
atically increase either ^^®Po or ^^^Rn counts. In this 
step, the fa has been measured, and the ion contribu¬ 
tion factor epo = 0.344 ± 0.003 is found using the Monte 
Carlo with a uniform initial z distribution of ^^^Rn. The 
expression relating these two activities is 

1214 

= (1 — /q + /aePo)(l — /a + /aepb)(l — fp + fpt-Bi) , 

^Rn 

where A^^/Ap^ = 0.116 ± 0.016. With this procedure 
fp is determined to be 76.3 ± 6.2%, which is in good 
agreement with the measurement starting from ^^®Po. 

Without considering the assumption that fa = fa, 
the primary source of uncertainty on fp (^^®Po —>.214 Pq^ 
is the ^^"^Po counting efficiency, contributing 2/3 of the 
total uncertainty. Uncertainties on epb and cbi are esti¬ 
mated by running the model with two extreme assump¬ 
tions for the unknown velocity of these ions, using vi and 
V 2 found in Section |IVB| as limits. The uncertainty on 
epo is calculated using the model by incorporating the 
uncertainties on the parameters G, ui, and V 2 - Finally, 
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treating the assumption that fa 7 ^ fay the limit that 
f'a = 0 the maximum value of fjs is 93%. Alternatively, 
for fp to be < 50% would require that > 76%. 


E. Alpha and beta ion and xenon hole production 

It is interesting to compare the alpha and beta decay 
daughter ion fractions to the surviving ion fraction of 
xenon holes, //j = 1 — r, where r is the fraction of initial 
ionizations that lead to recombination. For beta decay, 
recombination results in the generation of a photon, so 
the number of scintillation photons per ionization is r+a, 
where a is the initial ratio of the number of excited xenon 
atoms, or excitons, to ionizations. The ratio of photons 
[Nph) to electrons {Ne) from the decay is 

Nph _ a + r 
Ne ~ I-r 

which gives 


Using a mean experimental value for LXe, a = 0.13, 
pni^ and Nph = 30.5 keV~^ and = 42 keV~^ at the 
EXO-200 drift field from NEST [23] results of ^^"^Pb de¬ 
cay, fh = 0.65 for holes from beta decay. In comparison, 
for 5.3 MeV alphas, measurements give fh ~ 0.02 |M] . 
Since the ^^®Po daughter of ^^^Rn alpha decay should be 
located at the end of a 101 keV nuclear recoil track, it 
is perhaps more appropriate to compare to the hole ion 
fraction for a recoil track of such energy, fh = 0.18 pS] . 

The observed ion fraction for daughter ions is greater 
than that for Xe holes in both beta decay and alpha de¬ 
cay, as well as the nuclear recoil. This suggests the exis¬ 
tence of an additional reaction of the daughter atoms that 
competes with recombination, particularly in high den¬ 
sity ionization tracks in LXe. Both charge transfer from 
holes to the ^^®Po daughter atom and Penning ionization 
of the daughter atom by excitons (excited Xe atoms) are 
energetically allowed in liquid xenon and could substan¬ 
tially increase the daughter ion fraction relative to that 
of holes in alpha decay or nuclear recoil. These two pro¬ 
cesses could also contribute to enhancement of the daugh¬ 
ter ion fraction in beta decay, although to a lesser extent 
due to the lower density of ionization and excitation. 


V. CONCLUSION 

Alpha decays provide unique signatures in the EXO- 
200 detector that are useful to directly studying a number 
of low-energy physical processes. Detailed measurements 
of ^^^Rn and its decay products in the EXO-200 detector 
have been presented. 

Delayed coincidences between ^^^Rn and daughter 
^^®Po alpha decays have been used to determine the frac¬ 
tion of ^^®Po ions that result from ^^^Rn alpha decay in 
LXe, fa = 50.3 ± 3.0%, and to study the mobility of the 
daughter ion. Two distinct ion mobilities were found: 
0.390±0.006 cmV(kV s), and 0.219±0.004 cmV(kV s), 
with the ions initially moving at the higher mobility. The 
transition to the slower mobility was found to have a 
characteristic time that is proportional to the electron 
lifetime in LXe, suggesting that impurities play a role in 
the process. The fraction of ions from ^^"‘Pb de¬ 

cay in LXe was also determined to be fp = 76.4 ± 5.7%, 
assuming that the daughter ion fraction for the ^^®Po 
decay is the same as for the ^^^Rn decay. This result 
provides some basis for expecting high ionization frac¬ 
tion of ^^®Ba from double beta decay, which is relevant 
to the design of possible tagging methods. Progress on 
some of these tagging methods have been reported re¬ 
cently [a ID- The higher daughter ion fraction than hole 
ion fraction, particularly for alpha decay, was discussed, 
and a potential mechanism of charge transfer reactions 
with Xe holes and/or Penning ionization reactions with 
excitons was proposed to explain this. The EXO-200 ho¬ 
mogeneous LXe TPC has measured e~ drift velocities of 
1.7 X 10® mm/s [5], ion drift velocities of « 1 mm/s and 
neutral ion flows of « 10“^ mm/s, spanning a velocity 
range of 10®. 
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